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We explore the archetype problem of an escape dynamics occurring in a symmetric double well 
potential when the Brownian particle is driven by white Levy noise in a dynamical regime where 
inertial effects can safely be neglected. The behavior of escaping trajectories from one well to another 
is investigated by pointing to the special character that underpins the noise-induced discontinuity 
which is caused by the generalized Brownian paths that jump beyond the barrier location without 
actually hitting it. This fact implies that the boundary conditions for the mean first passage time 
(MFPT) are no longer determined by the well-known local boundary conditions that characterize the 
case with normal diffusion. By numerically implementing properly the set up boundary conditions, 
we investigate the survival probability and the average escape time as a function of the corresponding 
Levy white noise parameters. Depending on the value of the skewness /3 of the Levy noise, the escape 
can either become enhanced or suppressed: a negative asymmetry (3 causes typically a decrease 
for the escape rate while the rate itself depicts a non-monotonic behavior as a function of the 
stability index a which characterizes the jump length distribution of Levy noise, with a marked 
discontinuity occurring at a = 1. We find that the typical factor of "two" that characterizes for 
normal diffusion the ratio between the MFPT for well-bottom-to-well-bottom and well-bottom-to- 
barrier-top no longer holds true. For sufficiently high barriers the survival probabilities assume an 
exponential behavior. Distinct non-exponential deviations occur, however, for low barrier heights. 

PACS numbers: 05.10.Gg, 02.50.-r, 02.50.Ey, 82.20.Uv, 82.20.Wt 



I. INTRODUCTION 

The noise driven escape from a deterministically 
metastable state is a theme that impacts many phenom- 
ena in diverse fields of natural sciences Q, 0, El- H • In 
particular, as a notable model for a chemical reaction 
Kramers |l| pioneered the problem of an escape of a 
Brownian particle of mass m moving in a potential V(x) 
with local minima corresponding to an initial reactant 
and a final product state. In this scenario, both states 
are assumed to be separated by a barrier located at a 
position Xb- In the spatial-diffusion- limited regime, the 
Kramers rate theory is based on a stochastic dynamics 
that does not involve inertial effects and thus is described 
by an overdamped Langevin dynamics, reading 



dx 




(1) 
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Here, £(i) constitutes a white Gaussian noise process 
with correlations (C(i)C( s )) — 26(t—s), representing ther- 
mal fluctuations whose intensity is scaled by the friction 
•q. The escape problem then concerns the surmounting of 
an energetic barrier for stochastic trajectories that pre- 
dominantly dwell the neighborhood of separating attrac- 
tors which in this case are made of two neighboring po- 
tential wells. The imposed quasi-stationarity condition 
is realized by assuming special boundary conditions with 
respect to the time evolution equation for the probabil- 
ity density. In the overdamped regime, the evolution 
equation for the probability density p(x, t) follows the 
Smoluchowski dynamics 



dp(x, t) 
dt 



d_ 

dx 



J(x,t), 



where 



J(x,t) = — V'(x)p(x,t) 



ksT dp(x, t) 



V 



dx 
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and the stationarity approximation describes escape 
events which correspond to a constant, non- vanishing 
flux of probability J s — J(x). Those stochastic realiza- 
tions of the process that have surpassed the barrier top 
are immediately absorbed and re-inserted into the orig- 
inal attractor region. In this way a steady probability 
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flow across the activated barrier state located between 
the locally stable states of "reactants" and "products" 
is established. Escaped trajectories, absorbed at a posi- 
tion larger than the barrier location x a > Xb require that 
Ps{x&) = on the whole half- line x a > Xb- The rate for- 
mulation is then based on the "flux over the population" 
method 0, Q, yielding in this case of a Smoluchowski 
dynamics the celebrated result 0, 0, T 
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with AV denoting the activation energy AV = V{x\f) — 
V(x w ) and u) w , Ub representing the frequencies of the par- 
ticle's motion near the metastable potential minimum x w 
and at the top of the crossed barrier, respectively. This 
result represents a feasible estimate for the actual reac- 
tion rate, if all trajectories ejected by the source prop- 
erly thermalize before eventual thermal fluctuations drive 
them out of the initial well and, most importantly, a dis- 
tinct time-scale separation between escape dynamics and 
intra-well relaxation holds true. This latter requirement 
of a clear-cut time-scale separation is at the basis for the 
description of the escape dynamics in terms of a (time- 
independent) rate coefficient 0,0. 

Yet another alternative to the approaches discussed 
above is rooted in the concept of the mean first passage 
time (MFPT), i.e. the average time that a random walker 
starting out from a point Xq inside the initial domain 
of attraction, assumes in order to leave the attracting 
domain for the first time (2, tl L| ■ Put differently, the 
MFPT is the average time needed to cross the determin- 
istic separatrix- manifold for the first time @> EJ 0] • At 
weak noise the MFPT becomes essentially independent 
of the starting point, i.e. t(xo) ~ £mfpt for all starting 
configurations away from the immediate neighborhood 
of the separatrix. Given the fact that the crossing of the 
separatrix in either direction equals for normal diffusion 
one half, the total escape time t e equals 2£mfpt and thus 
the rate of escape k itself becomes in this case: 
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Characterization of the escape rate by use of the 
MFPT is a rather complex notion for a general class of 
stochastic processes. In particular, the MFPT analysis 
requires the choice of a correct boundary condition |a,l3- 
These are well known for one-dimensional stochastic dif- 
fusion Markov processes x(t) which are of the Fokker- 
Planck form, Eqs. (J2HHJ) or f° r one dimensional master 
equations with birth and death kinetics. With generally 
non-Gaussian white noise the knowledge of the boundary 
location alone typically cannot specify in full the corre- 



sponding boundary conditions for, say, absorption or re- 
flection, respectively 0,0- I n particular, the trajectories 
driven by non-Gaussian white noise depict discontinuous 
jumps. As a consequence, the location of the boundary 
itself is not hit by the majority of discontinuous sample 
trajectories. This implies that regimes beyond the loca- 
tion of the boundaries must be properly accounted for 
when setting up the boundary conditions. Most impor- 
tantly, returns (i.e. so termed re-crossings of the bound- 
ary location) from excursions beyond the specified state 
space back into this very finite interval where the process 
proceeds must be excluded. 

Following our reasoning in discussing Levy-Brownian 
motion on finite intervals |8( , we present in this paper the 
analysis of escape events of a non-inertial, generalized 
diffusion process which is driven by Levy noise dwelling 
a symmetric double- well potential. 



II. WHITE LEVY NOISE 

The standard definition of the Gaussian white noise 
specifies the latter as a time derivative of the Wiener 
process (i.e. a derivative of a stationary process with 
independent and Gaussian-distributed increments whose 
covariance is given by (W(t)W(s)) = min(t, s)). A classi- 
cal Brownian motion (the Wiener process) can be there- 
fore represented as a limit in distribution of independent 
Gaussian jumps taken at infinitesimally short time inter- 
vals of non-random length 1/n . Alternatively (following 
e.g. definition of Feller [ij 0, |nj), the Wiener process 
can be thought of as a limiting process of random Gaus- 
sian jumps at random Poissonian jump times: 



W(t) = UmYW h (±)= lim^f* 

n^oo ^ — ' \ n / n^oo y Ti 

N{tn) 

= lim V W k , (6) 

n — >oo ^ — » 



k=l 



The symbol [tn] stands for an integer number of jumps 
which for a Poisson counting process N(nt),t ^ with 
them mean (N(nt)) = nt are assumed to be independent 
(decoupled) of i.i.d random variables Wk sampled from 
the Gaussian distribution. The equality sign in Eq. © 
denotes a limit in distribution sense. For n — > oo, the 
Poisson distribution becomes peaked around k = nt and 
the limiting process W(t),t ^ tends to a Brownian 
motion diffusion for which the cumulative distribution 
function reads ^2] 
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Prob {W(t) < w} = Urn 5Z Prob (~ - w ) } P™b{N(nt) = k } 

/ w °° (nt) k 
V(27rfc/n)" 1/2 exp (-nx 2 /2k)^-— — exp (-nt)dx = 
■oot^c k\ 
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-dx. 
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Here, we introduce a non-Gaussian white noise as a 
derivative of the generalized Wiener process W a p(t), i.e. 
a non-Gaussian random process with stationary and in- 
dependent increments. The increments of such general- 
ized Wiener process have the a-stable distribution with 
the stability index a and the time increment At 1 / 01 as a 
scale parameter: 



t pt 
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where Q are independent random variables distributed 
with the stable, Levy probability density function (PDF) 
L a ,p{C\ a i A* = 0) an£ l NAs = t — tg. The parame- 
ter a denotes the stability index, yielding the asymp- 
totic power law for the jump length distribution being 
proportional to |C| — ■ The parameter a character- 
izes a scale, j3 defines an asymmetry (skewness) of the 
distribution, whereas fj, denotes the location parameter. 
Throughout the paper, we deal only with strictly stable 
distributions not exhibiting a drift regime; this implies 
a vanishing location parameter = in the remaining 
part of this work. For a ^ 1, the characteristic function 
(f)(k) — J x oo d(e~ lk( 'L a .i3((- 1 <j,fi) of an a-stable random 
variable £ can be represented by 



0(fc)=exp ~a a \k\ a (l - i/3sign(fc) 
while for a — 1 this expression reads 
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The three remaining parameters vary within the allowed 
regimes a 6 (0,2], j3 € [—1,1], a e (0, 00). Random 
variables £ corresponding to the characteristic functions 
© and (|TrHl can be generated using the Janicki-Weron 
algorithm [13| 

For a = 2 (and an arbitrary skewness parameter /3, cf. 
Eq. J3J), the generalized Levy-Brownian motion W a ^{t) 
becomes a standard Gaussian Wiener process whose time 
derivative leads to the well known Gaussian white noise 
limit. 

Following this interpretation, the stochastic source 
term in Eq. Q can be represented as a sum of indepen- 
dent pulses (having stable distribution) acting on equally 



spaced times. As a consequence, such Levy noise is white, 
i.e. its autocorrelation function is formally a Dirac-delta 
function of time. 

For the sake of clarification, we mention here that there 
exist different representations of other Levy noises that 
are non-Markovian in nature. As an example, so called 
fractional Gaussian (or Levy) noise is sometimes defined 
in literature as the time derivative of a fractional Brow- 
nian motion process 0, 0, 0, 0, 0] . In contrast to 
the Gaussian (or our Levy) white noise, the fractional 
Gaussian noise (fractional Levy stable motion) may ex- 
hibit slowly decaying time-correlations. This is however 
not the case with the type of Levy noise addressed in our 
paper, where the noise source at the level of the Langevin 
equation is a white noise process. 



III. ESCAPE IN A DOUBLE WELL: SURVIVAL 
PROBABILITY AND MEAN FIRST PASSAGE 
TIME 

Let us consider a Brownian, overdamped particle in an 
external potential V(x), i.e., 



V{x)=-\x* + \x\ 



(11) 



which is driven by Levy stable white noise. The sample 
trajectories are then obtained by a direct integration of 
Eq. JTJ: 



x(t)=x - / V'(x(s))ds+ / dL a .j3(s) 



l-o 



(12) 



to 



using the standard techniques of integration of stochas- 
tic differential e quat ion with respect to the Levy stable 
PDFs [1 [H Qll HHIi. The first passage time 
problem has been analyzed as r = inf{£ > | x(t) > Xb} 
with trajectories x(t) starting at x = xq and subject 
to corresponding boundary conditions. In particular, 
an absorbing boundary condition is realized by stop- 
ping the trajectory whenever it reaches the boundary, or, 
more typically, whenever it has jumped beyond that very 
boundary location. The role of reflection, which in the 
case of a free diffusion [g] has been assured by wrapping 
the hitting (or crossing) trajectory around the boundary 
location, while preserving its assigned length, see in Ref. 
0, |23 , is taken over naturally here by the confining po- 
tential walls of the symmetric double well. The details of 
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our employed numerical scheme for stochastic differential 
equations driven by Levy white noise has been detailed 
elsewhere [8|. In the following we shall omit cases when 
a = 1 with [3 0. In fact, this parameter set is known 
to induce instabilities in the numerical evaluation of cor- 
responding trajectories pi 111 El ElE^ . 

In the following we consider several differing situations. 
For a particle that starts out at a well bottom xq = —x w , 
cf. Fig. n and makes excursions toward the neighboring 
well bottom, taken as an absorbing boundary we use the 
subscript notation w — w. Likewise, fore a particle staring 
out at well bottom and being absorbed at barrier top 
location we use the notation w — b. Using the statistics 
of the firs passage time events we shall next study the 
two corresponding survival probabilities 

5 w _ b (<) = l-.F w _ b (t) (13) 
5 W _ W (<) = l-T w ^ w (t) , (14) 

where .F...(£) is the corresponding cumulative first pas- 
sage time distribution, i.e., 

/x b .x„ 
Pw— b,w— w{ x 7 t)dx , (15) 

with Pw-b.w-b(i) being the corresponding first passage 
time density. 

In order for the escape to become dominated by a clear- 
cut, single time-scale we use a sufficiently high potential 
barrier, thereby enforcing rare escape events. A too low 
barrier would involve too many recrossing events with 
the escape then being ruled by many time scales. Our 
used parameters for the potential are the parameter set 
a, b with a — 128,6 — 512, yielding a barrier height of 
AV = V(0) - V(-y^Jb) = = 8. This symmetric 
potential well is schematically depicted with Fig. ^ 
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FIG. 1: The generic double well potential V(x) = -%x 2 + jx 4 
for inspection of the Kramers problem discussed in the paper. 
To assure sufficiently high barrier separating the stationary 
states — x w ,a; w , the following set of parameters has been cho- 
sen: a = 128, b = 512, AV = V(0) - V^-y/aJb) = 8. 



A. Case of normal diffusion 

For a = 2, the process driven by white Levy sta- 
ble noise approaches the Brownian limit. The a-stable 
white Levy noise, as used in this study, then leads to 
the standard Gaussian white noise with intensity 2a 2 , i.e. 
(C(t)C( s ))a=2 = 2a 2 S(t — s). In this case, the correspond- 
ing values of the MFPT from the potential minimum — x w 
to the top of the absorbing potential barrier, i.e. Xb or 
to the neighboring, absorbing minimum, x = x w , can be 
calculated from the following quadrature formulas 

X 

MFPT(-x w -> x) = J cxp [V(z)/a 2 ] 

Z 

x J exp [-V(y)/cr 2 ] dydz. (16) 

— oo 

Note that for a system driven by white Gaussian noise 
the energy difference AV measured in units of ksT may 
be directly related to the intensity of the noise a 2 . In 
contrast, for Levy stable noises with a < 2, the scale 
parameter a is no longer "thermodynamically" related 
to the system temperature and consequently becomes a 
free parameter of the model. We have set throughout 
this study this value to a = y/2. 

B. Survival probabilities for a-stable noise driven 
escape 

Typical sample trajectories of the stochastic process 
defined by Eq. JIJ are depicted in Fig. [21 We observe 
that a non-zero (3 parameter induces a dynamical asym- 
metry for the escape dynamics occurring in a symmetric 
double well potential. It is reflected in our numerical sim- 
ulations by the emergence of stochastic trajectories that 
spend more time in the vicinity of one of the potential 
minima. As a consequence, one of the attracting states 
of the process (±£ w ) becomes favored over the other. 
This kind of behavior can be also detected in our other 
figures. Notably, for a decreasing value of the stability 
index a we observe larger fluctuations of the particle po- 
sitions. These occasional long jumps of trajectories may 
be of the order of, or even larger than the distance 2x w 
separating the two minima of the symmetric potential 
V(x). 

From the ensemble of single trajectories, which are sub- 
ject to the boundary conditions discussed above, and pre- 
sented in Fig. [21 we estimated the survival probability 
densities «S w _b(i) and £ w _ w (t), cf. see Fig. |3J The be- 
havior of the survival probability is consistent with the 
results as predicted by inspection of the corresponding 
stochastic trajectories: We clearly detect from Fig. [3] 
that at a chosen value of a a decrease of the skewness 
/3-parameter of a driving white noise causes a distinct 
decrease the rate of escape of the particle from the left 
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FIG. 2: Sample trajectories of random walks in the generic 
double well potential subjected to Levy stable noise with a = 
1.9 (left panel) and a = 0.9 (right panel) with various j3 (/3 = 
0, top panel, = —1.0, bottom panel). The scale parameter 
is a = v2, the chosen time step of integration is dt — 10 -4 . 



potential minimum; thus stabilizing the starting position 
at x = —x w . Most importantly, at sufficient high barrier 
heights a visible exponential decay of the survival prob- 
ability occurs, - according to S(t) = axp(— i/^MFPr)- 
This is the familiar behavior for regular Brownian mo- 
tion and has been detected already previously for sym- 
metric Levy noises [2J, |23 and for totally skewed (with 
(3=1) one sided Levy motions [26|. Here we observe 
it as well for skewed stable noises. Notably, the charac- 
teristic exponent Tmfpt clearly depends on both noise 
parameters, i.e. on the value of the stability index a and 
also on the skewness parameter (3. 

Our numerical analysis indicates that typically the sur- 
vival probability S w - W (t) for the "well-bottom-to- well- 
bottom" setup exceeds the survival probability S w ~b(t) 
for the "well-bottom-to-barrier-top" setup. This obser- 
vation can be explained by features of a noise-driven dy- 
namics: A particle performing the motion under the in- 
fluence of noise needs more time to dwell the neighboring 
potential minimum than is needed to reach the top of the 
barrier. This kind of behavior can be nevertheless weak- 
ened by diminishing the stability parameter a. For a < 2 
the trajectory of the random particle becomes discontin- 
uous and at sufficiently small a a particle, on average, 
completes its escape from the left potential minimum by 
a long jump which can overpass the right potential mini- 
mum. Consequently, both survival curves start to overlap 
converging to the same function at a given small stability 
index a. 

For the Cauchy limit, i.e., a = 1, an approximate re- 
sult for the mean crossing time as a function of the noise 
strength a has been given in Ref. [2j| by use of the 
fractional Fokker-Planck equation. In order to compare 
results obtained in this work with the former studies pre- 
sented by Chechkin et al. (24|, we have used a re-scaled 
form of Eq. Q 




20 40 60 80 
t 



20 40 60 80 
t 



dx 

17 



(x-x 3 ) + ((t), 



(17) 



FIG. 3: Survival probability densities S w -b(t) (left panel) 
and Sw-w(i) (right panel) for a particle dwelling the generic 
double well potential subjected to Levy stable noise with a = 
1.9 (top panel) and a = 0.9 (bottom panel). The simulation 
parameters are: dt = 10~ 5 , iV = 2 x 10 4 , a = \[%. 
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FIG. 4: Survival probability densities S w _b(t) and S w _ w (t) 
for a particle wandering in the generic double well potential 
subjected to Cauchy noise, i.e. Levy stable noise with a = 1.0. 
MFPT w _ b = 1.038 ± 0.008 and MFPT W _ W = 1.086 ± 0.008. 
The simulation parameters are: dt = 10 -5 , N = 2 x 10 4 , 
a = V2. 



obtained by a set of transformations x — > x/x w , t — > t/r 
with x^ — a/b, r — rj/a and consequently a a = D — > 
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t/x^. In these units a 



thus, 



n ■ , - w - 128 i/ 

an d On — ] = -^j. From the approximate formula derived 
in [24J for a strictly weak noise case at a = 1 one ob- 
tains MFPT pa 141.71 , whereas our simulations for this 
noise strength yield MFPT w 132.86, which differs only 
by ca. 7%. Exemplary survival probabilities S w -b(t) 
and S w - W (t) for this Cauchy-Levy noise are presented in 
Fig. 0] The data indicate that for the Cauchy noise the 
MFPT for well-bottom-to-barrier-top scenario is some- 
what smaller than the MFPT for the case of well-bottom- 
to-well-bottom, cf. Fig. 0] 

Our prior studies of a free Levy-Brownian motion 
driven by a stable noise has shown that for some pa- 
rameterizations of the Levy noise (totally skewed sta- 
ble distributions with a < 1) a non-exponential survival 
probability density emerges. Therefore, in this study we 
tested whether for a particle driven by a Levy-Smirnoff 
noise, i.e. with setting a — 0.5, = 1, a decrease of the 
potential barrier may induce visible deviations from the 
exponential distribution. The results of the test are dis- 
played in Fig. [S] where the survival probability densities 
S w -b(t) and S w ~ w (t) for various heights of the potential 
barrier AV are presented. For decreasing Ay the mo- 
tion of a particle indeed approaches the behavior of a free 
diffusion; consequently also distinct deviations from the 
exponential character of survival distributions do show 
up. 
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FIG. 5: Survival probabilities S'w-b(i) (left panel) and 
Sw-w(t) (right panel) for various barrier heights separating 
two wells of the potential. The system is driven by Levy- 
Smirnoff noise, i.e. Levy stable noise with a = 0.5 and /3 = 1. 



The simulation parameters are: 
a = V2. 



dt = 10"\ N = 2 x 10 4 



C. Behavior for the MFPT 

From the ensembles of collected first passage times 
we have also evaluated directly the mean values of the 
distributions. Our findings for MFPTs are presented in 
Figs. EMU The depicted results corroborate with the de- 
crease of the escape rate (i.e. the inverse of the MFPT) 



upon increasing the skewness (3. Moreover, for the stabil- 
ity index a = 2 (normal Brownian white noise) and for 
any skewness parameter (3, the Gaussian case is properly 
re-confirmed, cf. Figs.|HHZ| 





FIG. 6: Mean first passage times for well-bottom-to-barrier- 
top MFPT w _b (left panel) and mean first passage times 
for well-bottom-to-well-bottom MFPT W _ W (right panel) as a 
function of the skewness parameter /3. Simulation parameters 
dt = 1CT 5 , N = 2 x 10 4 , a = V2. The error bars estimated 
using bootstrap method with Nb = 2 X 10 3 . The white Levy 
stable noise with the stability index a = 2 and any allowed 
value of f3 is equivalent to the white Gaussian noise, what is 
manifested by the independence of MFPT(/3) for a = 2. The 
various symbols represent the various values of the stability 
index a: '+' a = 2.0, 'x' a = 1.9, 'o' a = 1.5, '□' a = 1.3, 
and '■' a = 1.1. 
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FIG. 7: Mean first passage times for well-bottom-to-barrier- 
top MFPT w _b (left panel) and mean first passage times 
for well-bottom-to-well-bottom MFPT W _ W (right panel) as 
a function of the stability index a. Simulation parameters 
dt = 10~ 5 , N = 2 x 10 4 , a = V2 and N b = 2 x 10 3 . The 
black squares indicate equality of MFPTs for a = 2 with any 
13. The various symbols represent the various values of the 
skewness parameter (3: '+' (3 = 1.0, 'x' f3 = 0.5, 'o' (3 = 0.0, 
and '□' (3 = -0.5. 
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D. MFPT-ratio no longer obeying normal 
behavior 

Next, we investigate the behavior of the value for 
the ratio R = MFPT w _ w /MFPT w _ b of the MFPT W _ W 
between the case with well-to well and well-to-bottom. 
Note that from the exact expression in Eq. (|16[1 . for a 
symmetric potential V(x) with normal Gaussian white 
noise fluctuations (i.e. for a = 2 and an arbitrary value 
of the skewness parameter (3) this ratio yields the com- 
monly known factor of R = 2. Put differently, in this 
case the random walker needs twice the time to reach 
the other potential minimum than to reach to the top 
of the potential barrier. In contrast, with a decreasing 
value of a the ratio R now distinctly deviates from 2, 
being always smaller than for a Gaussian diffusion case. 
Our numerical results are depicted in Fig. |SJ This de- 
viation can be understood by noticing that for a < 2, 
the stochastic escape trajectories of the random walks in 
the double well potential become discontinuous, meaning 
that the continuous movement of a particle becomes in- 
terrupted with long jumps. These occasional jumps are 
more probable to occur for small stability index a and ex- 
plain the observable statistics of the first passage times, 
implying in turn that R < 2. In particular, an increasing 
probability of long jumps over the barrier in the overall 
statistics of passages tends to equalize MFPT W _ W with 
MFPT w _b, yielding R values approaching 1. The effect 
of the skewness parameter (3 which is responsible for the 
asymmetric stochastic dynamics can be interpreted as 
an additional contribution to the directionality of ran- 
dom impact pulses which push the Levy-Brownian par- 
ticle within the potential well. For negative (positive) 
(3 w ±1 with a < 1, the driving Levy white noise be- 
comes a one-sided Levy process with strictly negative 
(positive) increments which stabilize trajectories around 
the starting position x = ±x w . 
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FIG. 8: Ratio R = MFPT w -w/MFPT w _ b of MFPTs to the 
neighboring potential well (MFPT W _ W ) and to the barrier 
top (MFPT w _b). The simulation parameters are: dt — 10 , 
TV = 2 x 10 4 , a = %/2 and N b = 2 x 10 3 . 



IV. CONCLUSIONS 

In this paper we have studied the survival probabilities 
and the mean values of the mean first passage times for 
escape from a symmetric double well potential when the 
overdamped dynamics is driven by general Levy white 
noise. The statistics of escaping trajectories is inves- 
tigated by a numerical analysis of a Levy noise driven 
Langevin equation with properly implemented boundary 
conditions. In contrast to former studies 0> HH aimed 
to understand statistics of barrier crossings with symmet- 
ric Levy flights or exploring the first passage problem for 
one-sided Levy motions 26], our work provides the first 
analysis of the Kramers problem with an arbitrary set 
of a, /3-parameters which characterize the white stable 
noise-source entering the Langevin dynamics. Following 
the Markovian character of the stochastic dynamics, - at 
sufficiently high barriers -, the time dependence of the 
survival probabilities within the potential well assume an 
exponential law. As can be expected, distinct deviations 
from the exponential behavior become detectable, how- 
ever, for low barrier heights when, similarly to a normal 
Gaussian diffusion case, many recrossing of the barrier 
are possible. 

An asymmetry of the Levy noise perturbing the dy- 
namics is shown to either enhance or also suppress the 
escape events and corroborates the intuitive influence of 
the (3 parameter which yields an additional, biasing con- 
tribution to the particles' motion due to the direction- 
ality of stochastic impact pulses. Moreover, the rate of 
escape is shown to exhibit a non-monotonic behavior as 
a function of the stability index a, with a discontinuity 
occurring at the value a — 1. 

In clear contrast to a normal diffusion behavior as 
typified by systems driven by white Gaussian noise, for 
which the random walker requires twice the time to 
reach the other potential minimum as compared to the 
mean time to reach to the (absorbing) top of the bar- 
rier, Levy white noise with a < 2 now causes a decrease 
in this ratio of the mean first passage times: The ratio 
R = MFPT w _ w /MFPT w _ b lies consistently below the 
"normal" value of R — 2. This remarkable result holds 
true for any chosen value of the skewness parameter (3. 

In summary, our numerical considerations demonstrate 
the richness of the bi-stable kinetics resulting from driven 
a-stable white noises. The observed features are rooted 
in the non-local jump lengths taken from a distribution 
which exhibit fat tails which in turn rule the random 
passages between the attracting states of a double well 
potential. 



Acknowledgments 

The Authors acknowledge the financial support from 
the Polish State Committee for Scientific Research 
(KBN) through the grants 1P03B06626 (2004-2005) and 
2P03B08225 (2003-2006), Marie Curie TOK COCOS 



8 



grant (6th EU Framework Programme under contract: 
MTKD-CT-2004-517186) and the ESF funds (E.G.N, 
and P.H.) via the STOCHDYN program. Addition- 
ally, B.D. acknowledges the financial support from the 
Foundation for Polish Science through the domestic 



grant for young scientists (2005). Computer simulations 
have been performed at the Academic Computer Center 
CYFRONET AGH, Krakow. P.H. also acknowledges the 
support by the German research foundation (DFG) via 
grant Ha 1517/26-1, 2 ("Single molecules"). 



[1] H. A. Kramers, Physica (Utrecht) 7, 284 (1940). 

[2] P. Hanggi, P. Talkner and M. Borkovec, Rev. Mod. Phys. 
62, 251 (1990). 

[3] P. Talkner and P. Hanggi, New Trends in Kramers Re- 
action Rate Theory, (Kluwer Academic Publ., Boston, 
1995); G. R. Fleming and P. Hanggi, Activated Bariier 
Crosing, (World Scientific, Singapore, 1993); P. Hanggi 
and J. Troe, Ber. Bunsen Phys. Chem. 95, 225 (1991). 

[4] P. Hanggi and P. Jung, Adv. Chem. Phys. 89, 239, 
(1995). 

[5] P. Reimann, G. J. Schmid and P. Hanggi, Phys. Rev. E 
60, Rl (1999). 

[6] C. Van den Broeck and P. Hanggi, Phys. Rev. A 30, 
2730 (1984); P. Hanggi, F. Marchesoni and P. Grigolini, 
Z. Phys. B 56, 333 (1984); P. Hanggi and P. Risebor- 
ough, Phys. Rev. A 27, 3379 (1983); V. Balakrishnan, 
C. Van den Broeck and P. Hanggi, Phys. Rev. A 38, 4213 
(1988). 

[7] P. Miiller, P. Talkner and P. Reimann, Physica A 247, 
338, (1997); P. Hanggi and P. Talkner, Phys. Rev. A 32, 
1934 (1985); 

[8] B. Dybiec, E. Gudowska-Nowak and P. Hanggi, Phys. 

Rev. E 73, 046104 (2006). 
[9] W. Feller, An introduction to Probability Theory and its 
Applications (John Wiley, New York, 1968). 

[10] V. M. Zolotarev, One Dimensional Stable Distributions 
(American Mathematical Society, Providence, 1986). 

[11] B. V. Gnedenko and A. N. Kolmogorov Limit Dis- 
tributions for Sums of Independent Random Variables, 
(Adison-Wesley, Reading, MA, 1954). 

[12] The result follows a general observation for separable con- 
tinuous time random walks (CTRW) which lead to nor- 
mal stochastic diffusion for finite mean waiting time for 



a jump to occur and finite step variance |lCj|. 

[13] R. Weron, Statist. Prob. Lett. 28, 165 (1996). 

[14] A. Janicki and A. Weron, Statistical Science 9, 109 
(1994); see also in: Simulation and Chaotic Behavior of 
a-Stable Stochastic Processes (Marcel Dekker, New York, 
1994); A. Janicki, Numerical and Statistical Approxi- 
mation of Stochastic Differential Equations with Non- 
Gaussian Measures, (Hugo Steinhaus Centre for Stochas- 
tic Methods, Wroclaw, 1996). 

[15] B. B. Mandelbrot, SIAM Rev., 10, 422 (1968). 

[16] S. Chaudhury and B. J. Cheravil, J. Chem. Phys. 125, 
024904 (2006). 

[17] M. Shlesinger, G. M. Zaslavsky and U. Frisch, Eds., Levy 
Flights and Related Topics in Physics (Springer Verlag, 
Berlin, 1994). 

[18] I. M. Sokolov, J. Klafter, Chaos, 15, 026103 (2005). 
[19] S. C. Kou and X. Sunney Xie, Phys. Rev. Lett. 93, 

180603 (2004). 
[20] P. D. Ditlevsen, Phys. Rev. E 60, 172 (1999). 
[21] A. Weron and R. Weron, Lecture Notes in Physics 457, 

379 (1995). 

[22] M. E. J. Newman and GT. Barkema, Monte Carlo 
Methods in Statistical Physics (Clarendon Press, Oxford, 
1999). 

[23] B. Dybiec and E. Gudowska-Nowak, Phys. Rev. E 69, 
016105 (2004). 

[24] A. V. Chechkin, V. Yu. Gonchar, J. Klafter and R. Met- 

zler, Europhys. Lett. 72, 348 (2005). 
[25] P. Imkeller and I. Pavlyukevich, J. Phys. A: Math. Gen. 

39, L237 (2006). 
[26] I. Eliazar and J. Klafter, Physica A 336, 219 (2004). 



